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Topological defects are ubiquitous from solid state physics to cosmology, 
where they drive phase transitions by proliferating as domain walls, monopoles 
or vortices. As quantum excitations, they often display fractional charge and 
anyonic statistics [Jj, making them relevant to topologically protected quantum 
computation [2^, but realizing a controlled physical resource for topological ex- 
citations has been difficult. Here we report evidence of topological excitations 
during the localization transition in strongly interacting two-dimensional elec- 
tron systems (2DESs) in GaAs/AlGaAs heterostructures. We find the electrical 
conductivity at low electron densities to follow a Berezinskii-Kosterlitz-Thouless 
(BKT)-like order-disorder transition implying a gate-tunable proliferation of 
charged topological defects [31, Hj. At low temperatures, a weakening in the 
temperature dependence of conductivity was observed, and linked to the zero 
point fiuctuations and delocalization of the defects. Our experiments also cast 
crucial insight on the nature of the ground state in strongly interacting 2DESs 
in presence of disorder 

At low temperatures (T), electrical conductivity a of a disordered 2DES decreases rapidly 
when the carrier density n^ is reduced. The dependence of a on Us often reflects the nature 
of the ground state which is electrically insulating. A classical percolation transition in 
inhomogenous 2DESs, for example, would cause a to vary as ~ {ug — n^Y^^ for Ug > uq, 
where Uq is the carrier density at the percolation threshold [3]. When disorder decreases, 
transition to the insulating state occurs at much lower rig, where the inter-particle Coulomb 
interaction is strong, and corresponds to a large interaction parameter r^ = l/a*^y/7Tn^ ^ 1 
{a*Q is the effective Bohr radius). The nature of the ground state in an interacting 2DES 
is far from clear, although possibilities of pinned charge-density waves (CDW), such as a 
Wigner crystal [6]-[T0]. charge stripes [HI [12], or other phases [13], have been suggested 
frequently. Physics of 2DES in this regime offers exotic possibilities such as metallic ground 
states, diverging effective mass and (^-factors, and so on [9]. 

In 2D, a CDW state with long range order does not exist at flnite T. However, states 
with quasi-long range order can exist, which can then undergo phase transitions into a 
disordered states through BKT type transitions involving the proliferation of topological 
defects [3], HI [TH - fTT] . Such transitions can in principle be effected by varying T or rig, and 
the correlation length (^) in the disordered state is a non-analytic function of the distance 



from the phase boundary in density and temperature space: 



^ oc exp 



(1) 



A 
_{X-XktY'^ 
where A is a constant, and X = n^ or T. When X > Xkt-, topological defects unbind, 

resulting in finite conductivity in the otherwise insulating pinned CDW state. The conduc- 
tivity varies in accordance with the number density riex (~ V'C^) of these defects, although 
an explicit evidence of this has not been observed until recently [18] . 

2DESs in modulation-doped GaAs/AlGaAs heterostructures are excellent testing ground 
for many-body physics, because it allows great control on the level of disorder. Similar to 
transverse magnetic field (-ffj.) fS], disorder can also stabilize CDW state, even in the absence 
of Hj^ , by suppressing the long wavelength fluctuations and arresting the kinetic energy of the 
electrons [lO]. This increases the melting point {i.e. tikt and Tkt), but excessive disorder 
may limit the order to very short range, or produce a glassy state, highlighting the need for a 
subtle balance. We use Si-doped GaAs/AlGaAs heterostructures, where disorder arises from 
(1) Coulomb potential of ionized Si-atoms, and (2) unintentional background doping. The 
former component was tuned by using devices from wafers with different setback distances 
{d) between the dopants and the hetero-interface (See Fig. la and Methods). It is important 
to prevent the device from pinching off {i.e. turn insulating) at high Ug which is the case in 
macroscopic systems with inhomogenous charge distribution fi9[ [20] . Mesoscopic 2DESs are 
more suitable since they are relatively insensitive to long range fluctuations in the conduction 
band [2n [22] , and hence form the backbone of our experiments. 

The devices were patterned in the form of a microbridge with a crossed metallic surface 
gate (Fig. la) (Ref [22]), limiting the effective area of the device to width of ~ 8 yum and 
length ~ 0.5 — 4 /im. With decreasing Us (negative potential on the gate), a dropped rapidly 
in all devices, as illustrated in Fig. lb and c. At low T (^ 0.5 K), small modulations on 
overall smooth variation of a were observed occasionally, which were not reproducible on 
thermal cycling. Moreover, a percolation transition also failed to describe the n^-dependence 
of a in our devices (dashed line in Fig. lb), indicating the influence of inhomogeneity to be 
small and sporadic. The pinch-off density where a — )■ 0, decreases with increasing T (shown 
in Fig. lb for device C2367), but increases with increasing H± (shown in Fig. Ic for T546 
and A2678). The density scale at which pinch-off occurs can be very low, n^ ^ 2 x 10^ cm~^, 
corresponding to r^ ~ 13, and hence the possibility of a CDW ground state arises naturally. 
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Figure 1 



FIG. 1: a) Schematic of a mesoscopic device geometry showing the 2DEG (red) at the hetero- 
interface and the dopant layer (light blue), b) Condutivity a (device C2367) as a function of 
electon density (n^) at different temperatures. The solid black lines are the fits for BKT transition 
(see Eq. 1 in text). The dashed line represents classical percolation description. The inset shows 
the linearity of {d\na /dns)~'^'^ in n^ {H_\_ = T, T = 3 K) confirming the BKT analysis. The 
x-intercept gives the critical density of melting {ukt)- c) Conductivity a of device T546 vs n^ at 
different magnetic fields and T = 262 niK. The fit illustrates the applicability of BKT description 
(Eq. 1 in text) over several decades in a. The inset shows similar fits for another device (A2678) 
with different disorder. 



A generic signature of pinned CDW, exploited for ID Peierls conductors [2^, as well as 
magnetically stabilized 2D electron solids [21], is the broad-band noise (BBN) noise from the 
electric-field induced thermally activated transitions between the metastable steady states. 
To explore this possibility, we have measured BBN in our mesoscopic devices in a 4-probe 
differential (ac + dc) mode, where the differential resistance of the mesoscopic region (dV/dl) 
was measured with a small low-frequency ac modulation (30 /iV) added onto a dc source- 
drain field {Vsd)- In the insulating regime (resistivity ^ h/e'^), dV/dl decreases, but without 
any abrupt threshold when \Vsd\ increases. This is illustrated for A2678 at T = 262 mK and 
H± = T in Fig. 2a. However, when dV/dl was recorded as a function of time (over several 
hours) at fixed Vsd, we observed a random switching noise for \Vsd\ ~ 100 — 200 /iV (see time 
traces in Fig. 2b). The net variance of noise [Sa'^/cr'^) hence shows a peak at a device specific 
Vsd of Vq ~ 150 — 200 /iV. Away from this regime, the noise is largely featureless with a 
power spectral density oc 1//" (Fig. 2c), with the spectral exponent a ~ 1 (Fig. 2d), arising 
from random potential fluctuations at the nearby acceptor sites. As shown in Fig. 2e, the 
peak in BBN decreases rapidly with increasing n^, and in the metallic regime, {a ^ e^//i), 
it decreases below the measurement background {6a'^/a'^ ^ 10^^°). The nonmonotonicity 
in noise is unlikely to arise from hopping, or electron glass behavior [25] , but can be taken as 
evidence of electrically driven dep inning transition in a pinned CDW [151 l2l]- The bias Vq 
at the noise maximum represents the depinning threshold which depends on the quenched 
disorder, and hence only weakly dependent of n^. The intermittent plastic flow of the CDW 
provides the switching noise for which a is expected to ~ 2 as indeed observed (Fig. 2d). 

Evidence of topological defects was subsequently obtained by focusing on a and scruti- 
nizing its dependence on Ug. In accordance with Eq. 1, we find a = ctq exp(— A/a/Ws — tikt) 
describes the variation of a over several decades in all devices over a wide range of T (base to 
4.5 K) and H± (see the fits in solid line in Fig. lb and c where Uo and A are constants). This 
can be readily associated with the proliferation of unbound charged topological defects which 
dominates electrical conduction in the disordered state (n^ > ukt)- The self-consistency 
of the fits was also confirmed by the linearity of {d\na /dns)^"^^^ in Ug (inset of Fig. lb), 
where the x-axis intercept gives the magnitude of tikt (within 10% of that obtained from 
nonlinear fit). 

To check if the analysis is physically meaningful, we constructed a melting phase dia- 
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Figure 2 



FIG. 2: a) Differential resistivity measurements (A2678) as a fuction of dc source-drain bias (Vsd) 
at 0.3 K and two carrier densities (10^ cm~^) corresponding to zero-bias resistivities of 4 /i/e^ and 
16 /i/e^. b) Time traces of normalized conductivity fluctuations at three different source-drain 
bias at Ug = 4.54 x 10^ cm^^. Data shows significant broad band noise at bias of Vq ~ 200 ^V. c) 
The power spectral density (PSD) of conductivity noise at the three bias showing a typical 1//" 
behaviour, where a is the spectral exponent, d) The plot of a as a function of dc bias which 
shows a peak magnitude of ~ 2 around Vq. The solid line is a guide to the eye. e) Normalized 
variance ba^ ja^ as a function of dc bias across the sample at different carrier densities (in units of 
10^ cm~^). The solid lines are guide to eye. The data shows a peak at around Vq signifying the 
depinning transition of the pinned electron solid. The peak value decreases at higher densities. 
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Figure 3 



FIG. 3: a) Critical melting density ukt (at two magnetic fields) for four samples with different 
setback distance d obtained from fitting the BKT (equation 1 text) to conductivity, ukt decreases 
with increasing d, except in C2367 (d = 60 nm) which has a strong background disorder (Methods). 
b) Evolution of rixT as a function of temperature at different magnetic fields {H±). c) Phase 
diagram of ukt in the T-H± space. The magnitude of ukt is indicated by the scale bar on 
the right. The quantum fluctuations, characterized by the parameter In /o-o which is the ratio of 
magnetic length to the inter-site distance, and the classical melting parameter F (see text) are 
evaluated at n^T and some of the traces are shown. A crossover from quantum to classical melting 
occurs at higher magnetic fields (bottom right) where F approaches ~ 80. 



gram with ukt extracted from the fits. Disorder quenches the kinetic energy of interacting 
electrons, and stabihzes the ordered state. Hence, decreasing the spacer thickness d gener- 
ally resulted in increasing n^r, unless the background doping dominates the disorder level. 
As shown in Fig. 3a, hkt decreases from A2407 {d = 20 nm) to T546 {d = 80 nm), but 
abruptly large for C2367 {d = 60 nm) most likely due to strong background doping (see 
Methods). In a given device though, disorder is fixed and ukt depends strongly on T and 
if_L, as illustrated for A2678 in Fig. 3b. Note that at H± = T, rixr is nearly 1 x 10^ cm~^ 
at low T, and corresponds to r^ ?^ 18 that is much smaller than the expected r^ for Wigner 
crystallization in 2DESs without disorder [B]. 

Fig. 3b illustrates that even a small H± can increase the melting point dramatically by 
quenching of quantum fluctuations of the constituent "atomic sites". It agrees with the 
recent suggestion that melting of a quantum electron solid needs to be addressed in terms of 



In/do rather than n^ alone [7]. (Here Ih = yh/eH±_ and Oq = 1/y^ are the magnetic length 
and mean electron distance, respectively). At higher if^ (~ 1.5 T), a further enhancement in 
riKT appears at low T, pointing towards a different mechanism governing the stability of the 
CDW. To understand this, we form a surface plot of ukt in T — H± space, where the traces 
of constant classical (F = e^ ^/rTs / 'iT[eoerkBT) and quantum {In/do) melting parameters at 
^s = Ukt are shown (Fig. 3c). The top-left corner, where both thermal and quantum 
fluctuations are strongest, the melting occurs at small ukt, while the bottom-right of the 
plot represents maximal stability. The abrupt increase in ukt for H± > 1.5 T occurs when 
Ih/o.0 ~ 0.3, signifying very little overlap of neighboring wavefunctions and a crossover 
from quantum to classical electron solid. Melting in this regime is described by F, which 
approaches ~ 80 in the extreme classical limit {H± > 2.5 T), agreeing closely with the 
estimation of F by Thouless for a Wigner crystal [14] . 

The T-dependence of a displays a peculiar weakening below a characteristic temperature 
scale O ~ 1 K (Fig. 4a-c) which needs further inspection. This behavior is common to 
all devices, and the magnitude of O was found to be insensitive to changing device length 
(over a factor of ~ 8), eliminating finite size effect as the cause. The dependence of a 
on Us also indicates Ug to be well-defined, and an explanation based on inhomogeneity- 
related models does not seem to apply here [HI [20]. To understand this quantitatively. 



we write cr{T) oc exp {—B/y/T — TxTjfiT). In a "Drude-like" scenario, for instance, 
cr(T) oc nex(T)T(T), where n^xiT) ~ exp [—Bj^JT — T^t) and /(T) would represent the 
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FIG. 4: a) Temperature dependence of conductivity a in three devices of different setback distances 
at various densities (in units of 10^*^ cm^^). The sohd hues are fits at the higher temperature range 
(T > 1.5 K) for temperature induced BKT transition (Eq. 1). d) ScaUng of scattering function 
/(T) obtained by dividing a by the BKT fit (see text) for different devices and densities. A 
single parameter scaling can be observed (holding the parameter for A2678 at 6.3 x 10^ cm~^ 
as reference). The line indicates an asymptotic behaviour of /(T) ~ 1/T^ for T <^ @. The top 
right inset shows experimental /(T) in various devices at different Ug (10^" cm^^). The bottom 
left inset shows the scaling parameter to increase with increasing disorder. 



temperature dependence of mean scattering time t{T). The behavior of f{T) changes at 
T ~ G, which was found to be close to the phonon gap of the electron solid through the 
following analysis: We obtain the domain size Lq {^ 430 nm) from the depinning threshold 
Vo (~ 200 /iV in BBN measurements for A2678 in Fig. 2) |2l], from which the phonon gap 
could be calculated as ~ hct/kBLo ^ 1.0 K (the sound velocity Ct ~ 4.7 x 10"^ ms~^ is taken 
to be that of a clean Wigner crystal) |26] . 

The behavior of cr(T) can be explained if T{T){r^ /(^)) is determined by an Andreev- 
Lifshitz (AL)-like mechanism of defect-mediated conduction in a quantum solid [2Z]. For 
T ^ 0, r — 7- constant, independent of T due to extensive phonon scattering, and the T- 
dependence of a is determined by that of nex{T) only. Hence to extract /(T), we fitted 
(t{T) by exp [—Bj^T — Tkt) in the range T > 1.5 K (lines in Fig. 4a-c) and subsequently 
divided it by the fitted trace. The increase in /(T) for T <^ 6 can be readily understood as 
freeze out of phonons and delocalization of defects as quantum objects with r ~ l/T^. A 
single parameter scaling of the T-axis collapses /(T) in all devices onto a universal curve, 
indicating B to be the only relevant energy scale (Fig. 4d). We find /(T) ~ r ~ 1/T^ at 
T ^ B (line) with B increasing with increasing disorder (i.e. decreasing d) (Fig. 4d(inset)). 

Summarizing, there are four key results of this paper: First, both noise measurements and 
scaling of a{n^ indicate a quasi ordered ground state in dilute 2DESs at low T. The precise 
nature of ordering is not known, Wigner crystallization being a strong possibility, although it 
is clear that disorder plays a crucial role in stabilizing such a phase. Second, the BKT scaling 
provides a new framework within which the transition from strong localization to metallic 
conduction can be viewed as an order-disorder melting transition. This is in contrast to 
percolation |5j, hopping [28j or Coulomb blockade [29] observed in many experiments over 
the years (mostly in larger or more disordered 2DESs). Third, the BKT transition also 
implies existence and proliferation of topological defects in 2DES in semiconductors for the 
first time with tunable number density. Finally, the scaling of the data described in Fig. 4d 
suggests that a mechanism like AL may apply here and that the topological defects behave 
as quantum mechanical entities capable of zero point motion and tunneling. 

METHODS: 
Wafer and device characteristics: 

The devices were fabricated from MBE (molecular beam epitaxy) grown 
GaAs/Alo.ssGao.eeAs heterostructures, where the two-dimensional electron gas was formed 
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300 nm below the surface. In all devices, delta (monolayer)-doping of Si (~ 2.5 x 10^^ cm~^) 
was implemented. The strength of the disorder is partially determined by the thick- 
ness of the spacer layer (undoped Alo.33Gao.66As) situated between the heterostructure 
interface and dopant layer. The typical as-grown electron mobility was found between 
0.9 - 3 X 10^ cmVV.s at n^ ^ 10^^ cm"^. The wafers T546 {d = 80 nm) and C2367 
{d = 60 nm) were grown in different chambers to observe the influence of background 
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0.5 


A2678 
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40 


1.8 


0.5 
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T546 
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80 
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Table 1. Details of devices used in our measurements showing the doping concen- 
tration fir, (10^^ cm~^), spacer thickness d (nm), as-grown mobility fi (10^ cm^/V.s) and 
background doping estimate Hbg (10^'^ cm~^). 
Measurements: 

The noise measurements were carried out inside a He-3 cryostat (base electron temperature 
Ri270 mK), using a four-probe ac-|-dc technique with the ac bias fixed at ^ 30 /iV at an 
excitation frequency of 128 Hz, and dc bias varied over ±400 fiV. The voltage fiuctuations 
across the sample were measured using a lock-in amplifier whose output was subsequently 
digitized using a 16-bit digitizer, followed by decimation and power spectral density (PSD) 
calculations using Welch's method of averaged periodogram. 
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